Spatial variation of fundamental couplings and Lunar Laser Ranging 
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If the fundamental constants of nature have a cosmic spatial variation, there will in general be 
extra forces with a preferred direction in space which violate the equivalence principle. We show 
that the millimeter-precision Apache Point Observatory Lunar Laser-ranging Operation provides a 
very sensitive probe of such variation that has the capability of detecting a cosmic gradient of the 
ratio between the quark masses and the strong interaction scale at the level V \n(m qusiT k/ Aqcd) ~ 
2.6 x 10 -6 Glyr" 1 , which is comparable to the cosmic gradients suggested by the recently reported 
measurements of Webb et al. We also point out the capability of presently planned improved 
equivalence principle tests, at the Ag/g < 10 -17 level, to probe similar cosmic gradients. 



1. INTRODUCTION 

Within many extensions of the Standard Model, the 
parameters of our fundamental theory need not be uni- 
versally constant but may vary in space and time. The 
search for such variations provides important constraints 
on such theories. Recently, Webb et al. [l| have reported 
evidence for a non-zero spatial variation of the fine struc- 
ture constant a. Parameterizing the variation of a by a 
dipole gradient 



a(x) 



(1) 



they find evidence, at the 4.2 a level, for a slope param- 
eter 



B a = (1.10 ±0.25) x 10~ 6 Glyr" 



(2) 



relative to the unit direction i a of right ascension a = 
17.4 ± 0.6 hours and declination 8 = —58 ± 6 degrees. In 
addition, Berengut et al. Q found weak indications for 
the existence of a gradient of the electron to proton mass 



ratio fj, 
slope 



m e /m f 



in the same direction z M = z c . 



Bf, = (2.6 ± 1.3) x 10" 6 Glyr" 



with 



(3) 



Other spatial gradients are much more weakly tested. 
For example, Donoghue and Donoghue Q have used the 
spatial constancy of the first acoustic peak in the Cosmic 
Microwave Background to bound a possible variation in 
the cosmological constant (or generalized dark energy) at 
the level of an analogous slope parameter 



B A < 0.91 x 10~ 2 Glyr~ 



(4) 



at the 95% confidence level. 

Because the masses of all the elements depend on the 
parameters of the Standard Model, a gradient in one of 



these parameters will lead to a force (as noted long ago 
by Dicke [![). Using the fine structure constant as an 
example, the dependence on a of the total mass-energy 
of system A, 



E A {a) = c 2 M A (a) , 



(5) 



implies that a spatial gradient Va of a , will lead to a 
force 



-VE A (a) 



2 dM A 
-c — — Va 



da 



(6) 



If we introduce the following dimensionless effective 
"charge" associated to the a dependence, 



Q a (A) = 



a 8M A 
Ma da 



(7) 



and parameterize the gradient of a by a slope and a unit 
direction as in Eq. ([T]), Va/a = B a z ai the above force 
reads 



-Q a (A) M A B, 



(8) 



If we now consider the dependence of the total mass- 
energy M A c 2 (in units of the Planck mass) of system A on 
the various dimensionless ratios (or coupling constants) 

7*2 — a, /i 5 ^quark 

/nip,. . . entering physics at energy scales 
< rnpC 2 , and if we assume the existence of (fractional) 
spatial gradients V In rj = B Ti z ri of the various dimen- 
sionless ratios, we see that body A will be submitted to 
an external acceleration, g A , of the form 

gA = w = ~ ? Qr * {A) Br * ° 2 2r * (9) 

i 

where Q ri = Q a ,Q^,--- are the various dimension- 
less effective "charges" associated to the dependence of 
the mass on the various ratios (or coupling constants), 
namely 



QrM) 



n 8M A d\n(M A /M P ) 



M A On 



d In Ti 



(10) 
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In the second form of the definition of Q ri we have re- 
called that Ma is to be expressed in units of the Planck 
mass Mp. [This corresponds to working in the "Ein- 
stein conformal frame" , where Newton's constant is held 
fixed.] 

If the various effective charges Q Ti {A) were indepen- 
dent of the considered body A, the result would be an 
unobservable (gravity-like) uniform free fall with a uni- 
versal acceleration g = g,4 = gs in a direction given 
by an average of the various gradients V In r j = B Ti z ri . 
However, composition dependence of (at least one of) 
the various charges, e.g. Q a {A) — Q a (B) ^ 0, or 
Q ll {A) — Q ll {B) ^ 0, will lead to differential accelerations 
gA — Sb 7^= and locally observable effects. Recently, we 
have studied @ the composition dependence of the ef- 
fective charges Q ri (A) corresponding to a complete set of 
dimensionless ratios entering low-energy physics, namely 
ro = Aqcd/Mp, and 

n = a,m u / A QC D,rn d / A QC D,rn e / Aqcd , (H) 

where m u , rrid are the masses of the light quarks. 1 Here, 
we separated the ratio ro — Aqcd/Mp, the depen- 
dence on which leads to composition-independent effects. 
We shall recall below our explicit results for the various 
charges Q ri (A) (i ^ 0). 

Each slope parameter B Ti defines a corresponding ac- 
celeration B ri c 2 , which enters the total acceleration 
multiplied by the corresponding dimensionless effective 
charge Q ri (A) . For instance, the a gradient ([2]) reported 
by Webb et al. [l| corresponds (using c/lyr= 950cm/s 2 ) 
to the acceleration level 

B a c 2 = (1.05 ± 0.24) x 10~ 12 cm/s 2 (12) 

while the acceleration level corresponding to the fi gra- 
dient suggested by Berengut et al. is 

B^c 2 = (2.5 ± 1.2) x 10" 12 cm/s 2 (13) 

The aim of this paper is to point out that the EP- 
violating effects of spatial gradients of a and of the 
mean quark-mass ratio 2 r^, = rh/ Aqcd (with rh = 
(m u +md)/2) at the levels Eqs. (fl"2"|) . (TT31 generate signals 
in the ranging to the Moon, which have a specific time 
structure, and an amplitude which seems large enough to 
be detectable by the recently started millimeter-precision 
Apache Point Observatory Lunar Laser-ranging Opera- 
tion (APOLLO) 0,0]. [See [1] for the results obtained 
from the pre- APOLLO Lunar Laser Ranging (LLR) ex- 
periments.] We also describe the weaker bounds on spa- 
tial gradients obtained by present laboratory-based ex- 



1 We have argued in that the composition dependence linked 
to the strange quark mass m a was subdominant. 

2 Note that, from a theoretical point of view, the ratio = 
m/h-QCD is akin to fi = m e /m p in the sense that it is the ratio 
between a lepton mass and an hadronic one. 



periments [lCf, and indicate that planned EP experi- 
ments at the Ag/g < 10~ 17 level will probe cosmic gra- 
dients at the levels of Eqs. JI2J), (fT5j) . 

Let us here emphasize the difference in outlook be- 
tween our previous work, and the present study. In 
Refs. H i|, we were considering the case where the 
spatial or temporal variation of a dimensionless parame- 
ter indicates the existence of a field, say ip, which carries 
the spacetime dependence, and we were considering the 
violations of the "weak version" of the Equivalence Prin- 
ciple (EP), i.e. composition-dependent accelerations of 
body A, mediated by the coupling of <p to local matter 
distributions. As a consequence, the locally observable 
EP-violating effects depended on the product of it two <p 
coupling strengths, say cxacie, where 

a A = d1n(M A ((p)/Mp)/dhLtp 

measures the coupling of ip to body A, and oie = 
d\n(ME{ t p)/Mp)/d\ntp its coupling to an external 
"source" body E (which could be the Earth, the Sun, 
or some laboratory source). However, we had normal- 
ized the definition of the fundamental couplings d ri of 
the "dilaton" field <p>, as they enter the low-energy La- 
grangian, so that we could write each a a in the specific 
form 

a A = d ro + ^2d n Q n (A) (14) 

exhibiting a simple factorization between the fundamen- 
tal dilaton couplings d ro , d Ti , and the phenomenological 
effective charges defined in Eq. (|10l) above. Note that 
there is no composition-dependent charge associated to 
the coupling to r = Aqcd/Mp, or, said differently, the 
charge Q ro {A) associated to r — Aqcd/Mp is simply 
Q ro (A) = 1 because, as the mass Ma can be written 
as the product of the hadronic mass scale Aqcd by a 
dimensionless function /(r^) of the dimensionless ratios, 
Eq. pip , the mass ratio Ma/Mp can be identically writ- 
ten as Ma/Mp = rof(ri). [Note also that the d Ti 's enter- 
ing Eq. (114(1 above correspond to the differences d Ti — d g 
in [5j, |6( (with d g = d ro ), because we defined above the 
ratios by Eq. pip which involved a logarithmic deriva- 
tive of Ma/ 'Aqcd, while we were working there with log- 
arithmic derivatives of Ma/Mp.] 

When contemplating, as we do here, possible varia- 
tions over cosmological distances, the field <p must be 
essentially massless. However, in the present work we 
shall not need to consider any specific model neither for 
the mass (or self-potential V((p)) of ip, nor for its mat- 
ter couplings d ro , d Ti . Indeed, the crucial point is that 
the observable acceleration © only depends on the ef- 
fective charges p"U|) , and on the various spatial gradient 
parameters V In = B Ti i ri . This makes the present 
investigation quite model independent, as well as inde- 
pendent from the usual interpretation of local EP tests 
(which involve the bilinear products (a a — cvb) &e)- 
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2. THE GRAVITATIONAL STARK EFFECT 
AND SPATIAL VARYING COUPLINGS 

The accurate monitoring of the lunar motion (most 
notably by LLR experiments |{|) has led to impressive 
tests of relativistic gravity, and notably of various as- 
pects of the EP |1JJ. Here we are interested in EP- 
violating effects in LLR that are linked to a fixed pre- 
ferred direction in space. Such effects have been studied 
by Damour and Schaefer (l2j in the context of binary 
pulsars. The analysis of such preferred-direction forces 
in the context of LLR has been done at leading order 
(LO) by Nordtvedt [Hf (see also [3]) , and to very high 
perturbative order by Damour and Vokrouhlicky [15| (us- 
ing the Hill-Brown lunar theory [.16]). For references to 
analytic studies of relativistic effects in lunar motion, as 
well as a self-contained introduction to Hill-Brown the- 
ory, see, e.g., [13] • Lunar dynamics is a notoriously dif- 
ficult problem because of the rather strong perturbation 
coming from the Sun's tidal forces, which leads to badly 
convergent perturbation series in powers of the parame- 
ter m = n'/(n - n') ~ 1/12.3687. [Here, n' = 2tt/(1 yr) 
denotes the mean sidereal angular velocity of the Earth 
around the Sun, and n — 27r/(27.32 days) [l8| the mean 
sidereal angular velocity of the Moon around the Earth.] 
For some effects, a LO perturbation treatment in m can 
be significantly inaccurate both because of the occurrence 
of small denominators, and of the slow convergence of 
the m-perturbation series. This is, for instance, the case 
for the Laplace-Nordtvedt effect of polarization of the 
Moon's orbit by an EP-violation linked to the Sun's grav- 
ity where higher-order terms in m increase the leading- 
order result by more than 62% [U El- 

In the case of 

interest here of what has been called the "gravitational 
Stark effect" [l2[ , i.e. the perturbing influence of a differ- 
ential force (with a fixed direction) acting on a gravita- 
tionally bound two-body system, the situation is similar, 
though with significant differences. 

Let us first recall that the classical (electric or gravita- 
tional) Stark effect is an example of singular perturbation 
where a small perturbing force can have a large effect. 
If we were approximating the dynamics of the relative 
Earth-Moon coordinate r = xm — xe in the presence 
of an external acceleration Ag = gj\/ — gE by means 
of the Lagrangian [after factorization of the Earth-Moon 
reduced mass [i= itie tum /{piE + itim)] 



1 f dr\ G(m E +m M ) , . 
L= 2U + r + A9 ' r ' (15) 



we could find the exact solution of the perturbed dy- 
namics by separating the Hamilton- Jacobi equation cor- 
responding to the Lagrangian (|15[) in parabolic coordi- 
nates (£ = r + z, rj = r — z, </>), with a z axis oriented 
along Ag. One then finds that the exact solution corre- 
sponding to elliptic orbits undergoes a complicated sec- 
ular evolution during which the osculating elements of 
the elliptic motion wander very far away from any given 



initial state. For instance, even if the perturbing accel- 
eration Ag is very small, the osculating eccentricity will 
not undergo small oscillations around its initial value eo 
but will, on time scales na/Ag, take values quite differ- 
ent from eo- This instability of elliptic motion under a 
constant force can also be seen by using the averaged evo- 
lution equations of the osculating orbital elements. More 
precisely, if we consider the evolution of the semi-major 
axis a, of the Lagrange-Laplace-Runge-Lenz eccentricity 
vector e — ea (directed towards the periastron), where 
(a, b, c) are orthornormal unit vectors with a pointing 
towards the periastron and c along the orbital angular 
momentum £ = (1 — e 2 ) 1 / 2 c = r x v, one finds averaged 
evolution equations of the form [l2| 



da 
y dt' 

K dt' 
M 

K dt' 



, 

f x £ , 
f x e 



where 



f = 3Ag 
2 na 



(16) 



(17) 



with n denoting the sidereal angular frequency of the 
Moon. 

We see that while a stays secularly constant, the vec- 
tors e and £ rotate one into another. More precisely (with 
/ = /z), one easily sees that, while e 2 and i z stay con- 



stant, the two complex combinations e a 



rotate as 



- „+*/* . 



s (0); e v (t) 



-ift 

P J ; 



i L, and 



,(0), (18) 



e^ ~\~ £y 



and |e„| 2 = 



leaving constant \e z 
(consistently with e 2 + £ 2 = 1 = const.). 

This Stark instability is rooted in the well known de- 
generacy of the Coulomb problem, i.e. the fact that the 
radial uj t and angular frequencies happen to be ex- 
actly equal, uj r — lo^ = n, for a 1/r interaction potential. 
As a consequence, any lifting of the Coulomb degener- 
acy by an additional interaction potential (causing tu r to 
differ from to^) will tame the Stark instability. Ref. [l2[ 
considered the case where this lifting was due to the gen- 
eral relativistic modifications of the 1/r Newtonian po- 
tential. In that case the only modification of the secular 
evolution equations Eqs. (|16l) is the appearance of an 
additional contribution +td p c x e on the r.h.s. of the 
evolution equation of e, where oj p = oj^ — ui r is the pre- 
cession frequency of the binary system, due to relativistic 
effects. 

As a first orientation towards understanding the Stark 
effect in the lunar motion, let us start by assuming that, 
as in the case studied in Ref. |12j , it is enough to replace 
the second secular evolution equation in Eqs. (|16|) by 



,de 



f x £ + oj v c x e 



(19) 
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where w p describes the precession of the orbit of the 
Moon, which occurs with a period of 8.85 years [l8[. [Let 
us note in passing that the amplitude of the eccentricity 
evolves according to 4| = (1 — e 2 ) 1 ^ 2 f • b.] 

In addition, as the eccentricity of the Earth- Moon sys- 
tem is small e = 0.0549, we can, in first approximation 
(in view of the appearance of e on the r.h.s. of the evo- 
lution equation of the angular momentum £) , neglect the 
small wobbling of the direction of the orbital angular 
momentum c and approximate equation (| 19[) (using also 
£ = (1 - e 2 ) 1 /^ ~ c) by 



de 

( — ) = f x c + cjpCX e = Up c x (e — e/) 



where we have introduced 



3 A gJ 



Up 2 naujp 



(20) 



(21) 



where f j_ = f — (f • c) c is the component of the external 
force in the plane of the orbit. 

The general solution of Eq. (|2T)]) reads 



(22) 



where the constant vector ej describes a "forced eccen- 
tricity" (or a "polarization") induced by the external 
EP-violating force, and where e p (t) is the usual, u p - 
precessing free eccentricity, which is allowed in absence of 
external Stark effect. Note in passing that the polariza- 
tion of the elliptic orbit by the external force is oriented 
in the opposite direction of the projection fj_ cx Agj_ of 
the force on the orbital plane. Indeed the eccentricity 
vector defines the direction towards the periastron, so 
that the elliptic orbit is mainly elongated in the apoas- 
tron direction — ey cx — fj_. 

The final observable result, of this leading order treat- 
ment, linked to the polarization e/, is a sidereal frequency 
oscillation of the Earth-Moon range, connected with the 
direction of the projection fj_ of the perturbing accelera- 
tion onto the orbital plane, of the form 

A LO r{t)=p^°cos{n(t-t Q )-<j> f ) . (23) 

where the (algebraic) amplitude is 



f±a _ 3 Ag ± 
u n 2 nuj r , 



(24) 



and where n(t — to) is the longitude of the Moon, and 
4>f is the longitude of the direction fj_/|fj_|, both lon- 
gitudes being measured within the orbital plane, from 
some common origin. [Neglecting the small inclination 
i ~ 5 degrees of the Moon's orbit on the ecliptic, we 
can consider that both longitudes are ecliptic longitudes, 
counted from the vernal equinox.] 

The above result was obtained as a leading-order ap- 
proximation, under the simplifying assumption that the 
main effect of the solar tide on the Earth-Moon system 



was to introduce a perigee precession term in Eq. (|19l) . 
However, the solar tide has more effects than this, and, as 
we recalled above, the theory of the lunar motion under 
the combined effect of the G(uie + rriM)/r Earth-Moon 
potential and of the quadrupolar tide ^x l x^ dij(Gms / D) 
of the Sun, is a notoriously subtle theory, exhibiting many 
instances of slowly converging perturbation expansions in 
the (not so small) expansion parameter 



0.080849 : 



12.3687 



(25) 



Although the basic small parameter entering lunar the- 
ory is the ratio of the solar tidal potential to the Earth- 
Moon potential, which is of order m 2 ~ 10 -2 , the lu- 
nar perturbation theory is not an expansion in powers 
of m 2 ~ 10~ 2 , but proceeds (beyond the LO term) in 
powers of m ~ 1/12 (see, e.g. [16[)- We note, in particu- 
lar, that the perigee precession frequency of the Moon is 
given by a perturbation expansion of the form 



w. 



177 

IF' 



1659 



85205 



3073531 
" 2 13 -3 



which converges so slowly that the sum of higher-order 
terms approximately doubles the LO analytical result 
ujp° /n — 3m 2 / 4 (see [2(| for the literal computation of 
the perturbation expansion of uj p /ti up to the eleventh 
power of m). Let us also note that the above LO re- 
sult for the range perturbation due to the Stark effect 
contained Lu p /n = 3m 2 /4 + ■ • • as a small denominator 
that significantly amplified the effect of the external per- 
turbing force f . [The presence of a small denominator is 
linked to the instability, recalled above, of elliptic motion 
under a constant force, because this small denominator 
tends to zero in the limiting case of the Lagrangian (TT5"j) .] 

The appearance of such small denominators oblige 
one to tackle in a more complete manner the Stark ef- 
fect on the Moon's orbital motion. This was done by 
Damour and Vokrouhlicky [l| using a Hill-Brown treat- 
ment (with the help of a dedicated algebraic computer 
programme). More specifically, Ref. [15| worked out the 
lunar range perturbation Ar(t) induced by an external 
acceleration Ag to very high order in the powers of m. 
This range perturbation Ar(t) is the sum of many dif- 
ferent frequency components that come from the non- 
linear combination of the basic sidereal frequency n of 
the Moon (linked to the angular distance between the 
Moon and the external fixed direction Ag, or rather 
its projection Ag ± on the Moon's orbital plane), with 

linked 



even multiples of the synodic frequency n - 
to the angular distance (seen from the Earth) between 
the Moon and the Sun. Among the spectrum of com- 
bined frequencies ±(n + 2j(n — n')) (j £ Z), two of them 
were found to be dominant: the basic sidereal frequency 
in (of period 27.32 days), and the j = — 1 combination 
±(n - 2(n - n')) = t(« - 2n') (of period 32.13 days). 
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The result of [l~5j can be written 3 as 



Ar(t) = Pf 



cos [n(t — to) — (j>f] 



15 S'Am) 
■ m — r cos [n(t — to) — 2r(t) — <j>t\ 



8 S f (m) 



"26) 



Here rit) = (n — n')t + to denotes the synodic phase, i.e. 
the angular distance between the Moon and the Sun, and 
the overall amplitude is given by 



(27) 



where A.gj_ and <pf are the magnitude, and the longitude, 
of the projection Agj_ of the external acceleration onto 
the lunar orbital plane , and where Sf(m) and S'Jm) 
are m-perturbation series that start as 1 + O(m). For 
instance, the beginning of the expansion of St(m) reads 



Sf(m) 



1 - ^ m + ^fp m 2 



127 637 3 
384 " l 



4 172 299 4 
2304 



0(m 5 ) 



(28) 



and Table IV of [15| gives the coefficients of this expan- 
sion to the ninth order in m. Even with such a high- 
order expansion one finds that the last term is still of 
fractional order 1CP 3 . This slow convergence is related 
to the presence of a pole in the series Sf(m) and S'Am) 
near m cr ~ —0.18407. To get an accurate numerical esti- 
mate of these series Ref. [15j used a Pade resummation, 
with the results [for m — tomooii given by Eq. (|25jl ] 



S f (m) ~ 0.5050 



and 



15 S' f (m) 
— m 



S f (m) 5.94 

for the fractional coefficient of the subleading term with 
frequency n — 2(n — n') = — (n — 2n'). 

An observationally important aspect of the result (j2"6")l 
is the appearance of a specific combination of two har- 
monics, with known periods and phase, and with nearly 
comparable magnitudes. In particular, the fact that the 
amplitude of the n — 2n' harmonic is only ~ 6 times 
smaller than the LO n harmonic is a result of the sub- 
tleties of lunar perturbation theory. This term comes 
from the basic solar tide perturbation which is propor- 
tional to m 2 , but it has been amplified to the 0{m) 



3 Here we change the notation of Section IV of |15| : e.g. 
A G N G Ag ± , SgaiM ^ S f (m), etc. 

4 Even in the full Hill-Brown treatment (with neglect of the lu- 
nar inclination, and of the eccentricity of the Moon's orbit), one 
finds that only the projection Agj^ of the external acceleration 
matters for the range perturbation. 



level by a small denominator (with the additional factor 
15/8 ~ 2, leading to 15to/8 ~ 1/6). We have checked the 
presence of this subleading term by directly solving the 
forced Hill's equation [MEIIIoi d 2 q(T)/dT 2 +Q(T)q(T) = 
o~(t) for the transverse perturbation g(r) to Hill's varia- 
tional orbit. Here, r = (n— n')t+To as above. In that for- 
mulation, the frequency component qj exp[i(l + m + 2j)r] 
comes with the denominator 0q — (1 + m + 2j) 2 which is 
small when j = — 1, being 0(m) [namely, Qq — (m— l) 2 = 
[1 + 2m + 0(m 2 )] - (m - l) 2 = 4m + 0(m 2 )]. 

For what concerns the leading term, with frequency 
n, in Eq. (|2T)]) . it corresponds to the result Eq. (|2"5)l 
of the approximate treatment explained above. In both 
cases (f>f is the longitude of the external acceleration pro- 
jected within the orbital plane. Note that if one were 
using the LO analytical results {uj p /n) LO = (3/4)m 2 and 
Sf°(m) = 1, Eq. ([24]) would read -2Ag±/n' 2 and would 
agree with Eq. (|27p. However, the exact result Eq. (|2"?| 
is smaller than this by about a factor two, because of the 
correcting factor Sf(m) ~ 0.5050. As noted in Ref. (l5j . 
when using in Eq. ([M]) the actual perigee precession 
uip (which is about twice larger than its LO estimate 
(uj p /n) LO — (3/4)to 2 ), one captures most of the effect 
of the slowly converging series Sf(m). 

We can finally apply the result Eq. (|2"oT) to the EP- 
violating acceleration Ag = gM — Se where each is 
given by equation ©. Let us first note that the result 
only depends on the amplitude and longitude of the pro- 
jection on the orbital plane of the vectorial sum of the 
external accelerations 



Ag = g M - SE = - $3(Qr«(M) " Qr<(E)) B r 



. (29) 

Alternatively, one could write the range perturbation as 
a sum of terms of the type of the r.h.s. of Eq. (|26l) . 
each one having an amplitude Pfivi) and a phase <j)f(ri). 
[Note that we consider here algebraic amplitudes (that 
can be negative), and that the longitudes 4>f{ r i) always 
refer to the direction of the projection i ri ^_ of the gra- 
dient direction +z ri .] Let us, as it simplifies the writing 
of our results, make the natural assumption that all the 
cosmological gradients of the coupling constants are par- 
allel to each other, i.e. z ri = z independently of the label 
ri, and let us denote by <pf the longitude of the projected 
gradient direction zj_. This leads to a total range per- 
turbation of the form Eq. (f2"6")> . with a total (algebraic) 
amplitude of the form (after cancellation of two minus 
signs) 

pf = 2S f (m) £(Q n (M) - Q ri (E)) (30) 



where B Tt ± is the magnitude of the projected gradient 
B r . zj_ . 

Note that the numerical prefactor 2Sf(m) ~ 1.010 
is close to 1, and that the parameter combination 
B ri ±c 2 /n' 2 , where we recall that n' — 27r/lyr, is the 
product of an acceleration by the square of a time, and 
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is indeed a length. [Alternatively, we can think of it as 
the product of the spatial gradient B riJ _ — [length]" 1 , by 
the squared length c 2 /n' 2 = (1 lyr) 2 /47r 2 ]. As a first ori- 
entation, note that a gradient of order B ~ 10~ 6 Glyr - , 
i.e. Be 2 ~ 10 -12 cm/s (similar to the recent suggestions, 
Eqs. ©, 0, dT2j), (HSJ)) corresponds to a figure of merit 
Be 2 /n' 2 ~ 25 cm. In order to estimate the corresponding 
signal in lunar laser ranging, we need next to estimate the 
numerical value of the various EP- violating charge differ- 
ences Q ri {M) — Q ri {E) corresponding to the difference 
in composition of the Moon and the Earth. This will be 
the focus of the next Section. 

Before tackling this issue, let us briefly mention that 
the central values of the equatorial coordinates a = 261 
degrees, S — —58 degrees of the cosmological gradient 
of the fine structure constant reported in (lj corresponds 
to an ecliptic latitude equal to (3 — —34.7 degrees, and 
an ecliptic longitude equal to A = —95.8 degrees. The 
latter ecliptic longitude predicts 5 the value of the longi- 
tude entering the range perturbation Eq. (f2l)|). namely 
4>f = A. On the other hand, the ecliptic latitude /3 enters 
the observable range pf through the projection of the 
cosmological gradient direction onto the orbital plane, 
i.e. (essentially) onto the ecliptic. More precisely we have 
B ri ± = B Ti cos/3. Note that cos/3 = 0.822 for the gradi- 
ent reported by Webb et al. so that this projection (that 
we shall include in the estimates of the next Section) re- 
duces only by 18% the full possible observable effect of 
such a gradient on the lunar motion. 



eration ([9]) will come from the dependence of Ma/ Aqcd 
on the masses of the light quarks, on the electron mass 
and on the electromagnetic fine structure constant. For 
the quark masses, instead of considering separately the 
masses of the up and down quarks, it is convenient to 
consider their average and difference, namely 



m 



(m d + m u )/2, dm = (m d - m u ) (31) 



and to work with the fine structure constant a together 
with the dimensionless ratios 



A 



QCD 



A 



QCD 



A, 



QCD 



(32) 



Because the fermion masses are the product of a dimen- 
sionless Yukawa coupling Tj and the Higgs vacuum ex- 
pectation value (vev) v, rrii — Tiv/\/2, these ratios can be 
considered as the product of the dimensionless I\ by the 
ratio v/ Aqcd between the basic weak-interaction scale v 
and the basic strong- interaction scale Aqcd- In view of 
the independence of the mechanisms leading to the ap- 
pearance of the two basic scales v and Aqc d , it seems a 
priori theoretically natural to expect that the cosmolog- 
ical gradients (if any) of the three mass ratios (|32p will 
be of similar magnitudes, and therefore similar to that 
of the ratio /i = m e /m p for which the value ([3]) has been 
recently suggested. 

Refs. [5J, |6| derived the folllowing approximate esti- 
mates for the four effective charges Q ri associated to the 
three mass ratios (I32D. and to a: 



3. EP VIOLATING CHARGES 

Let us now estimate the numerical values of the various 
EP- violating charge differences Q ri (M) — Q ri (E) enter- 
ing the magnitude of the cosmologically induced Earth- 
Moon differential acceleration @. This issue has been 
discussed in our previous work [3, H| where we described 
the leading dependence of atomic masses on the parame- 
ters of the Standard Model. Because of the large neutron 
and proton masses, the dominant determining factor for 
all atomic masses is simply the scale of the strong interac- 
tions Aqcd- However, as already mentioned above, the 
main dependence of atomic masses on Aqcd is universal, 
and is conveniently factored out by rewriting each atomic 
mass Ma as Ma = AQC£>/(rA, where the fj's (for i ^ 0) 
denote the dimensionless ratios Eq. (fTTj) . As a conse- 
quence, the charge Q ro associated to r$ = Aqqd/Mp via 
the general definition Eq. (TTU|) is simply Q ro (A) = 1, and 
the only non-zero contributions to the differential accel- 



5 As above, we are here neglecting the small inclination of the 
Moon's orbit on the ecliptic. Note that taking into account this 
inclination, and its secular variation, would predict an additional 
small adiabatic variation of the amplitude pf (and phase (pf) 
of the range perturbation Ar(t), corresponding to the secular 
wobbling of the lunar orbital plane. 



Qrh 



0.093- ^?-0.020-^ 2Z ' 2 



1.4 x 10 



A 1 / 3 

-iZ(Z-l) 



A 2 



Qsm — Fa 



Qm c — Fa 



0.0017 ■ 



A 4 /3 



A-2Z 



A 



5.5 x 10~ 4 - 
A 



and 

Q a = F A 



-1.4- 



4 



7.7 



Z(Z-l) 



A 4 / 3 



x 10" 



(33) 



(34) 



(35) 



(36) 



where Fa = Am amu /MA, with m amu — 931 MeV denot- 
ing the atomic mass unit. The common factor Fa is very 
close to 1, and we shall replace it by 1 in our estimates 
below. 

Approximating the Moon as made of silicate (i.e. es- 
sentially SiC>2), and the Earth as made of a mantle of sili- 
cate and a core of iron (representing 32% of its mass), the 
above formulas yield the following Moon-Earth charge 
differences 
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Qrh{M)-Q A (E) 

Q8m(M) - Qsm(E) 
QrnA M ) ~ QrnA E ) 

Q a (M)-Q a (E) 



-1.1 x i(r 3 

-3.8 x 10~ 5 

1-6.1 x icr 6 
-3.i x icr 4 



(37) 
(38) 
(39) 
(40) 



We see that the most important charges for EP vio- 
lation are Q m and Q a . This dominance of Q„ and Q a 
over the other charges is true for most values of Z, A. 
It is related both to the rather small coefficients enter- 
ing Qsm, Eq. (j34f . and Q mc , Eq. (|35|) . and to the fact 
that A ~ 2Z along the periodic table. If we use the 
facts that F A - 1 = O(10~ 3 ), and that A ~ 2Z, we can 
simplify the composition dependence of the main EP- 
violating charges Q m and Q a as Q m ~ Q^+cst. and 
Qa — Q' a +cst., where [1, H| 



Qm 

Q'a 



0.036 



E4 \- 10-* (41) 

-7.7^7^ x 10- 4 (42) 

We list the values of — Q' m and Q' a for a sample of ma- 
terials in Table I. This list shows that the maximum 
value of a charge difference would be the Q m differ- 
ence between a heavy element and a light one with 
(heavy) — Q A (light) ~ + 10 -2 . The corresponding 
maximum acceleration level AgJ? ax = AQ m ax B m c 2 ~ 
10~ 2 B m c 2 would numerically be 



10- 6 Glyr" 1 



x 10" i4 cm/s 



for a cosmic gradient of order of those reported by Webb 
et al. 



TABLE I: Approximate EP-violating 'effective charges' for 
a sample of materials. These charges are averaged over the 
(isotopic or chemical, for SiC>2) composition. 



laterial 


A 


Z 


Qrh 




( 






Li 


7 


3 


18.88 xlO 


-3 


0.345 


xlO 


-3 


Be 


9 


4 


17.40 xlO 


-3 


0.494 


xlO" 


-3 


Al 


27 


13 


12.27 xlO 


-3 


1.48 


xlO" 


-3 


Si 


28.1 


14 


12.1 xlO" 


3 


1.64 


xlO" 


-3 


Si0 2 






13.39 xlO 


-3 


1.34 


xlO" 


-3 


Ti 


47.9 


22 


10.28 xlO" 


-3 


2.04 


xlO" 


-3 


Fe 


56 


26 


9.83 xlO" 


3 


2.34 


xlO" 


-3 


Cu 


63.6 


29 


9.47 xlO" 


3 


2.46 


xlO" 


-3 


Cs 


133 


55 


7.67 xlO" 


3 


3.37 


xlO"" 


-3 


Ft 


195.1 


78 


6.95 xlO" 


3 


4.09 


xlO" 


-3 



4. OBSERVATIONAL SIGNALS LINKED TO 
POSSIBLE COSMOLOGICAL GRADIENTS 



reported values 6 Eqs. ©, © [and (13]), (fl2])]. we 
conclude that those cosmological gradients entail EP- 
violating differential accelerations on the Moon directed 
along the unit vector zj_ (with ecliptic longitude equal to 
(j)f ~ A ~ —95.8 degrees), and with algebraic magnitudes 



(43) 



where 



Ag m _ 



and 



= +2.2 



-{Q m {M) - Q m {E)) B m± c 2 
B m 



(44) 



2.6 x 10- 6 Glyr~ 



— x 10 _i;, cm/s 



-(Q a (M) - Qa(E)) B a± c 2 
-2.7 ' 



(45) 



1.1 x lfr 6 Glyr^ 1 



x l(r lb cm/ 



These differential accelerations then induce a perturba- 
tion in the Earth- Moon range which has the specific time 
signature (I26|) with an algebraic magnitude 



where 



and 



Pa = 



Pf = Prh + Pa 



2.6 x 10- 6 Glyr" 



Ba 

1.1 x 10- 6 Glyr" 



0.59 mm 



0.068 mm 



(46) 



(47) 



(48) 



As we see, the dominant effect is expected to be linked 
to the cosmological gradient of iti/Aqcd, so that the re- 
cent findings of Webb and collaborators [H, 0] suggest 
the presence of millimeter-level sidereal fluctuations in 
the Earth-Moon range. Such fluctuations, if they ex- 
ist, should be detectable by the APOLLO experiment. 
Indeed, this experiment has shown its capability of ob- 
taining "normal-point" range measurements with nightly 
median uncertainty of 1.8 mm for their entire data set, 
and 1.1 mm for their recent data [8j. By accumulat- 
ing millimeter-level range data over a sufficiently long 
time (comparable to the perigee period ~ 8.85 yr), the 
APOLLO experiment should be able both to decorre- 
late the specific sidereal signal (|2"6"|) from the many ex- 
isting Newtonian range effects (which include synodic, 
n — n', and anomalistic, w r = n — ui p , frequencies), and 
to measure its amplitude pf to a fraction of a millime- 
ter. Depending on the result of such an analysis, the 



Putting together our results, keeping only the domi- 
nant terms linked to Q m and Q ai and scaling the possi- 
ble cosmological gradients of fh and a by the recently 



As mentioned above, it is natural to expect that a cosmo- 
logical gradient of fj, = m e /m p implies a similar gradient in 
the weak-interaction/strong-interaction ratios Ti, and notably in 
m/A Q cD- 
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APOLLO experiment could either establish the reality 
of a cosmological gradient of coupling constants, or set 
upper bounds on the gradients and B a (or more pre- 
cisely on the combination B^i + 0.28B Q j_ entering Ag±) 
at levels smaller than the levels Eqs. ([3]), ([2]) [and (fl"3|). 

Let us mention in this respect that a sidereal range 
perturbation of the approximate form (|23[) . with a phase 
<f> f linked to the center of the Galaxy, has been searched 
for in the pre-APOLLO, few-centimeter-level LLR data 
after the suggestion of Ref. 13:]. Nordtvedt, Miiller and 
Soffel [5l[ published an upper limit of 



Affgal < 3 x 10 



-14 / 2 

cm/s 



(49) 



on a possible perturbing differential acceleration linked 
to the Galactic center, while a further analysis of Miiller 
(cited in (22|) obtained 



A# ga i = (4 ±4) x 10~ 14 cm/s 



(50) 



The gain in sensitivity of the APOLLO experiment, by 
more than a factor 10, and the richer time signature of the 
more complete signal (|26p , make us expect that it will be 
possible to probe the acceleration levels (j4"4")l . (|4"5")l . above. 

It is instructive to compare the (potential) sensitivity 
of LLR experiments to external EP-violating accelera- 
tions with the sensitivity of other EP experiments. There 
have been constraints on anomalous cosmic accelerations 
by laboratory based EP tests. The most precise is quoted 
as a differential acceleration in any direction of the sky 

m 

\g(Be) - g(Ti)\ < 8.8 x 10" 13 cm/s 2 (95%C.L.) 

The charge differences Q Ti (Be) — Q Ti (Ti) are again domi- 
nated by Qrh(Be)-Qrn(Ti) = -7.23x 10~ 3 and Q a (Be)- 
Qo,(Ti) = —1.56 x 10~ 3 . Assuming, as above, that the ef- 
fect of the gradient of to (which couples to the dominant 
charge difference) dominates, the above upper bound on 
|<?(Be) — <?(Ti)| can be readily converted to a bound on 
the corresponding cosmological gradient Bm , namely 



\B fn \ < 1.3 x 10~ 4 Glyr" 



(51) 



This is is weaker than the recently suggested (theoreti- 
cally similar) gradient Eq. (|3|) by a factor ~ 50. Such a 
difference in acceleration sensitivity between Earth-based 
EP experiments and LLR ones might seem surprising 
in view of the fact that both types of experiments cur- 
rently lead to comparable limits on the (Eotvos) EP- 
violation parameter n = Ag/g, namely ?7EarthMoon = 
(-1.0±1.4)xl0- 13 H, versus r? B eTi = (0.3± 1.8) x 10~ 13 
[lfj | , and that both types of experiments use comparable 
background accelerations g in the ratio Ag/g. [Indeed, 
the g due to the Sun at Earth is gs — 0.6 cm/s 2 , while 
torsion balance experiments use only the horizontal com- 
ponent of the Earth gravity, namely gs _i_ — 1-7 cm/s 2 at 



a latitude of 45 degrees.] We note that the greater sen- 
sitivity of LLR experiments to external (especially fixed- 
direction) accelerations is essentially rooted in the spe- 
cific Stark instability mentioned above. Indeed, generally 
speaking, a differential acceleration Ag acting during a 
characteristic time t c (which is t c ~ lu^ 1 = T/(2w) for a 
periodic phenomenon of angular frequency us and period 
T) corresponds to a measurable displacement of order 
Ar ~ Agt 2 — Ag/u 2 . In the LLR case, we saw above 
that the range perturbation is Ar ~ Ag/n' 2 which is 
larger than the expected perturbation ~ Ag/n 2 associ- 
ated to the lunar frequency n by a factor 



lyear 



lsidereal month 



(13. 37) 2 = 178.7 



This amplification factor lies at the root of the in- 
creased sensitivity of LLR experiments to external EP- 
violating accelerations having a fixed direction. We note 
in passing that the LLR sensitivity to the usually con- 
sidered Laplace-Nordtvedt solar-rooted EP-violating ac- 
celeration is only amplified, w.r.t. Ag/n 2 , by a paramet- 
rically smaller factor ~ (3/2) (n/n') ~ 20, i.e. about ten 
times less than in the "Stark" , fixed direction case. This 
difference is due to the difference in the corresponding 
nearly resonant denominators, namely, in Hill's equation, 
a denominator 6> — (1 + m) 2 = [1 + 2m + 0{m 2 )] — (1 + 
to) 2 = 0(to 2 ) in the sidereal-frequency (Stark) case, ver- 
sus 6» - (l) 2 = [1 + 2to + 0(to 2 )] - 1 = 2to + 0(m 2 ) in 
the synodic-frequency (Laplace-Nordtvedt) one. 

Let us finally note that presently planned improved EP 
tests such as the Satellite Test of the Equivalence Prin- 
ciple (STEP) [H (r) ~ 10" 18 ) or proposed cold-atom- 
technology tests ((n ~ 10 -17 ) (24j, which will make use 
of the full strength of the Earth gravity, gE — 980 cm/s 2 , 
will be able to probe the cosmological-gradient-induced 
differential accelerations discussed above. Indeed, the ac- 
celeration (|4"4")1 linked to a cosmic gradient of to can be 
rewritten (modulo a cosine factor, with a specific time 
dependence linked to the projection onto the sensitive 
direction of the considered EP experiment) as 



n 



— = 2 5 9E ( A< ^™ 



xlO 



-17 



g g V 10~ 2 J \ 2.6 x lO" 6 Glyr" 

where we allowed the considered man-made EP test to 
optimize the choice of materials by having a AQm ~ 
10~ 2 , i.e. ten times better than for the Earth- Moon case 
(see Table I). This result shows that the LLR test of the 
cosmological acceleration (|4"4")) . that should be doable by 
the APOLLO experiment, corresponds (from the point of 
view of the sensitivity to a cosmic gradient) to an Earth- 
based EP test at the rj ~ 10~ 17 level. 
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